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Abstract. We present a temperature map and a temperature profile of the central part (r < 20' or 1/4 virial radius) 
of the Coma cluster. We combined 5 overlapping pointings made with XMM/EPIC/MOS and extracted spectra 
in boxes of 3.5' x 3.5'. The temperature distribution around the two central galaxies is remarkably homogeneous 
(r < 10'), contrary to previous ASCA results, suggesting that the core is actually in a relaxed state. At larger 
distance from the cluster center we do see evidence for recent matter accretion. We confirm the cool area in the 
direction of NGC 4921, probably due to gas stripped from an infalling group. We find indications of a hot front 
in the South West, in the direction of NGC4839, probably due to an adiabatic compression. 
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1. Introduction 

Numerical simulations of cluster evolution (e.g. Evrard 



199C; Schindler & Muller 1993) suggest that the temper- 



ature structure of the Infra-Cluster Medium (ICM) is a 
powerful indicator of the evolutionary state of clusters. 
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In particular the accretion of a sub-cluster, a common 
phenomenon in standard hierarchical formation scenario, 
should manifest itself by characteristic features in the tem- 
perature map, like heated gas between the two units just 
before the collision. 



Recent studies of the Coma cluster with the ASCA 
1999| , Watanabe et al. |1999[ ) 



satellite (Donnelly et al. 
revealed complex temperature variations in this mas- 
sive cluster. They were interpreted as indicative of re- 
cent mergers, confirming earlier evidence based on optical 
dynamical studies (Colless & Dunn 1996 and references 
there in) and X-r ay m orphological an alysis (Bri el et al. 
1992t White et al. |l993| , Vikhlinin et al. |l994| , |l997|) . ASCA 
covered a broad energy band, which is essential for pre- 
cise temperature estimate, but the observations suffered 
from a relatively large energy dependent PSF. Therefore 
temperature structure determination with ASCA might 
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Table 1. Observations 



Obs. 


Rev 


RA. 


DEC. 


MOS1 


MOS2 




MOS1&2 counts 








(J2000.0) 


(J2000.0) 


Exp. 


Exp. 


[0.3-10] 


keV 


[5-10] 


keV 










(ksec) 


(ksec) 


Source 


Bkgd 


Source 


Bkgd 


Pc 


86 


12 h 59 m 47 s 


27°57'00" 


16.4 


16.3 


7.70 10 b 


4.2 10 4 


3.52 10 4 


1.00 10 4 


P5 


86 


12 ft 59 m 28 s 


27°46'53" 


20.9 


21.4 


6.63 10 5 


5.5 10 4 


3.23 10 4 


1.30 10 4 


P6 


93 


12 h 58 m 50 s 


27°58'52" 


7.4 


7.3 


1.81 10 5 


1.9 10 4 


8.53 10 3 


4.52 10 3 


P9 


93 


13 ft 00 m 33 s 


27°56'59" 


20.8 


21.0 


7.03 10 5 


5.4 10 4 


3.08 10 4 


1.29 10 4 


P10 


98 


12 h 59 m 38 s 


28°07'40" 


20.6 


20.9 


5.26 10 5 


5.4 10 4 


2.43 10 4 


1.27 10 4 



have been subject to systematic errors. Furthermore the 
spatial resolution was insufficient to resolve precisely the 
temperature radial profile in the very core of clusters. 



The EPIC instrument (Turner et al. 2001) on board 
XMM (Jansen et al. |200l| ) combines a high sensitivity 
with good spatial and spectral resolution, on a wide energy 
range. In this paper, we use this unique capability to study 
the temperature structure in the central (8 < 20' = 0.78 
Mpc) region of Coma. We present further XMM results 
in two other papers of this issue: the large scale morphol- 
ogy of Coma (Briel et al. |2001| ) and the dy namic s of the 
infalling NGC 4839 group (Neumann et al. |200l| ). In the 
following, we assume Hq = 50km/s/Mpc and qo = 0.5. 



2. Data Analysis 

2.1. The data 

The central part of Coma was observed with 5 overlap- 
ping pointings in Full Frame mode with the EPIC/MOS 
camera (medium filter) and in extended Full Frame mode 
with the pn camera. As CTE correction in this pn mode 
is still being studied, we considered only MOS data in the 
present spectroscopic analysis. 

We generated calibrated event files with SASv4.1, ex- 
cept for the gain correction. Correct PI channels are ob- 
tained by interpolating gain values obtained from the 
nearest observations of the on board calibration source. 
Data were also checked to remove any remaining bright 
pixels. We excluded periods of high background induced 
by solar flare protons. We discarded all frames correspond- 
ing to a count rate greater than 15 ct/lOOs in the [10 — 12] 
keVband, where the emission is dominated by the particle 
induced background. Finally, only events in the nominal 
FOV are considered. 

The central position of each pointed observation is 
listed in Table 1, together with the revolution number and 
remaining observing time after cleaning. 

2.2. Spatially resolved spectroscopy 

Spectra in various regions were considered to study tem- 
perature variations. Each pointing and each MOS camera 
are first treated separately. 



2.2.1. Correction for vignetting effects 

The effective area at a given energy depends on position. 
When extracting the spectrum of a region, we weight each 
photon falling at position {xj,yj) of the detector and of 
energy Ej by the ratio of the effective area at that posi- 
tion, to the central effective area (for this energy). This 
weighting is taken into account in the error estimate. The 
'corrected' spectrum obtained is an estimate of the spec- 
trum one would get if the detector was flat. A detailed 
description of the method and of the vignetting calibra- 
tion data used are given in Arnaud et al. ( [2001 ). Note that 
the vignetting due to the RCA is included, but is assumed 
to be energy independent (the variations are less than 1% 
below 6 keV). 



2.2.2. Background estimate 

We generated EPIC MOS background event files (one for 
each MOS camera) by combining several high galactic 
latitude pointings. The data are cleaned for bright pix- 
els, background flares and regions corresponding to bright 
point sources are excluded. The integrated exposure time 
is 94.3 ksec for MOS1 and 78.9 ksec for MOS2. These event 
files can be used for a proper estimate of the cosmic ray 
(CR) induced background but not the X-ray background, 
which depends on pointing position and filter. However, 
bright cluster emission, like the one observed in the Coma 
center, usually dominates the backg round except at high 
energy (see also Arnaud et al. 2001 ). We thus are mostly 
sensitive to CR induced background. Note also that the 
offset pointings considered here always include the cluster 
peak emission, so scattered light is not a problem. The to- 
tal estimated number of source and background photons, 
in the [0.3 — 10] keV and [5 — 10] keV energy ranges, are 
listed in Tab.|l] for each pointing. Note the degradation of 
the S /N ratio at high energies as a consequence of the very 
hard CR induced background. 

It is known that the CR background changes slightly 
in the FOV. It is thus better to consider the same extrac- 
tion regions in detector coordinates for the source and the 
background. Furthermore, if one wants to combine spectra 
of a given physical Coma region obtained from different 
pointings, it is mandatory to define extraction regions in 
sky coordinates. To alleviate this practical problem, we 
simply generated a specific background event file for each 
Coma pointing and camera by modifying the sky coordi- 



Arnaud et al.: The temperature structure in the central region of Coma 



3 




Fig. 1. The EPIC/MOS mosaic image of the central region of 
Coma (5 overlapping pointings) in the [0.3 — 2] energy band. 
The iso contours are the residuals (in a) after subtracting the 
best fit 2-D f3 model. The step size is 4cr and the lowest iso- 
contour corresponds to 3a significance. The position of the 
bright galaxies are marked. 



nates of the background event file using the aspect solution 
of the considered Coma observation. 

For consistency, the background spectra were obtained 
using the same correction method as used for the source. 
The background component induced by CR is not vi- 
gnetted, but as we extract the background and source 
spectra from the same region in detector coordinates the 
correction factor is the same and does not introduces bias. 

Further details on the characteristic of the EPIC/MOS 
background and subtraction method can be found in 
Arnaud et al. 



2001 



2.2.3. Spectra extraction and spectral fitting 

The source and background spectra of a given region (de- 
fined in sky coordinates) are first extracted separately for 
each MOS camera and each pointing, using the method de- 
scribed above. As the spectra are corrected for vignetting, 
the spectra of the same physical region observed with dif- 
ferent off-axis angle (from different pointings) and differ- 
ent camera can be simply added to maximize the signal 
to noise ratio. The errors are propagated using quadratic 
summation. 

Before model fitting, the source spectra are binned so 
that the S/N ratio is greater than 3 a in each bin af- 
ter background subtraction. The spectra are fitted with 
XSPEC using isothermal mekal models (with fixed red- 
shift z = 0.0231). Although the thermodynamic state of 
the plasma could be more com plex ( e.g. see the isobaric 
multiphase model of Nagai et al. 2000 ), this is adequate to 
study temperature spatial variations, the derived best fit 
temperature being an estimate of the mean temperature 




"0.9 r 



channel energy (keV) 

Fig. 2. EPIC/MOS1 (green) and EPIC/MOS2 (blue) spectra 
extracted from within 10' in radius of the galaxy NGC 4874. 
Red line: best fit redshifted isothermal model: kT = 8.25 keV 
and an abundance of 0.25. Bottom panel: residuals between 
model and data. 

in each projected region considered. Since the spectra are 
'corrected' for vignetting effects we can use the on axis 
MOS response file, which is considered to be the same for 
MOS1 and MOS2 (version v3.15). Only data above 0.3 
keVare considered due to remaining uncertainties in the 
MOS detector response below this energy. Unless other- 
wise stated errors are with a 90% confidence level. 



2.3. Imaging analysis 

The MOS mosaic image in the [0.3 — 2] keV energy band is 
presented in Fig. [l]. The count images for each camera and 
each pointing are extracted using the weighting procedure 
described above to correct for vignetting. They are then 
projected on a common sky reference axis, summed and 
divided by the mosaic exposure map (which takes into 
account the exposure time and FOV coverage for each 
pointing). The images are not background subtracted. In 
that energy range and in this central area of the cluster, 
the particle background is negligible for MOS. 

3. Results 

3.1. Core morphology 

To identify significant substructures in the core, we fitted 
the MOS image with a 2-D ellipsoidal (3 model plus back- 
ground and built up the map of the residuals of the data 
over the best fit model. The method is discussed in detail 
in Neumann & Bohringer ( 1997| ) and Neumann (1999). 
The iso-contours of significance of the excess (number of 
er over background plus cluster model) are overlaid on the 
MOS image in Fig. |. 

We unambiguously confirm the excess emission around 
the two central galaxies (NGC4874 and NGC 4889) and 
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Fig. 3. Position of the 3.5 x 3.5' boxes, in which we derive 
EPIC/MOS spectra. A different symbol is used for each of the 
5 pointed observations. When the same region is observed in 
two or more pointings, the corresponding spectra are summed 
(see text). The spectrum of each region is fitted with an isother- 
mal model to build the temperature map displayed in Fig. [E|. 
Circles: central 9 < 10' FOV of each pointed observation. 

the tail of gas in the direction of NGC4911, revealed by the 
wavelet analysis of Vikhlinin et al. (1997). However, con- 
trary to Vikhlinin et al. result, this hlamcntary structure 
does not seem to be directly connected to the Coma center, 
but originates 0.5 Mpc South of it. Note also that part of 
the excess is due to resolved galaxy emission (NGC 4911, 
NGC 4921, QS01259+281, and 5 other point sources). 

Diffuse excess emission is clearly detected in the South- 
West in particular in the direction of the NGC4839 group 
(see Briel et al. 2001, for full discussion). 
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Fig. 4. Temperature determined from the offset pointings (P5, 
P6, P9 and/or P10) versus the temperature obtained, for the 
same region, using the central pointing (Pc) data. Errors are 
at la-level. The insert shows the histogram of the differences 
in term of \ 2 ( see text for details). 

21cm value and x 2 is unchanged. In the following we thus 
fix the ./Vh value to the 21 cm value. 

When fitted separately the MOS1 and MOS2 tempera- 
tures are consistent and the relative normalization is 1.05. 
In the following analysis we will thus sum MOS1 and 
MOS2 spectra. 

The best fit overall temperature in this central (R < 
10') region is in remarkable agreement with the overall 
GINGA value of kT = 8.21 ± 0.09 keV (Hughes et al. 



1993) and is marginally consistent with the ASCA value 



of 9. ± 0.6 keV ( Donnelly et al. |1999[ ) obtained for the 
central (R < 9') region. Note that the hard excess seen 
by SAX (Fusco-Femiano et al. 1999) could not be seen by 
XMM, since it appears above 20 keV. 



3.2. Overall spectrum 



3.3. Temperature map 



To compare with results obtained with other satellites we 
extracted the overall MOS1 and MOS2 spectra from a cir- 
cular region of 10 arcmin in radius centered on NGC4874. 
Only data from the central pointing are used. The data 
are fitted in the [0.3 — 10] keV range, with independent 
normalizations for MOS1 and MOS2, and common tem- 
peratures and abundances. The spectra are plotted on 
Fig. |2[ the bottom panel gives the residuals. Fixing the 
(7V H value to the 21 cm value (iV H = 8.95 x 10 19 cm -2 , 
Dickey & Lockman 1990| ) , we obtain a best fit temperature 
of kT = 8.25 ± 0.10 keV and an abundance of 0.25 ± 0.02. 
The reduced X 2 is 1-38 ( X 2 = 1457 for 1058 d.o.f). The fit 
is satisfactory. The residuals are concentrated around the 
instrument edges, with residual ratios between data and 
model of about ±5%, consistent with our present knowl- 
edge of the instrument response (Fig. ||). If we let the A^h 
value free we get kT = 8.20 keV, an abundance of 0.25 
and A^h = 9.4 ± 0.9 x 10 19 cm~ 2 , in agreement with the 



To study the cluster temperature structure, we next ex- 
tracted spectra in 3.5' x 3.5' contiguous regions in sky 
coordinates. The box size was chosen so that the two cen- 
tral galaxies fall approximately in the center of a box, 
and that a sufficient S/N ratio is reached for each box. 
Circular regions (20" in radius) around bright sources (in 
particular NGC 4889 and NGC 4911) are excised from 
the boxes. The overall region considered for this spatially 
resolved spectroscopic analysis is about 20' in radius. We 
only considered boxes at off-axis angles smaller than 10' in 
each pointing (the vignetting factor being more uncertain 
beyond this radius). The central 6 < 10' region of each 
pointing, delineated as circle, and the central position of 
the various boxes (95 in total) are plotted on Fig. ||. 

The validity of our vignetting correction can be as- 
sessed by comparing the fitted temperature of the same 
region in various pointings. The vignetting effect (decrease 
of effective area with off-axis angle) increases with en- 
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Fig. 5. Color coded temperature map. Note the hot front in 
the south-west (white) and the cold region in the South-East 
(brown/dark red). The isocontours of the PN image in the 
[0.3 — 2.] keV band (Briel et al. 2001) are superimposed. The 
lowest contour corresponds to 6.3 x 10 -3 ct/s/arcmin 2 and the 
step size is 4. x 10 -3 ct/s/arcmin 2 



ergy. An understimate (overestimate) of this energy de- 
pendence would yield to underestimate (overestimate) of 
the temperature. Since the various pointings of the same 
region correspond to different off- axis angles, an improper 
vignetting correction would translate in systematic differ- 
ences between temperature estimates for the same region. 
Fig. |J shows the temperatures determined from the offset 
pointings T Q ff versus the temperature obtained, for the 
same region, using the central pointing T c (errors are at 
l<7-level) . The insert shows the histogram of the differences 
in term of X 2 = (T oS - T c ) 2 /{a(T oS f + <r{T c ) 2 ) computed 
for each pair of measurements. Three outliers (\ 2 > 4 or 
more than 2a discrepancy in estimates) are clearly appar- 
ent. They correspond to the points at (6. 7, 8. 4), (6. 5, 7. 3) 
and (7.6,8.5). We do not see any particular problem in 
the corresponding spectra (the statistic and fit are good) 
and failed to find any obvious reason for the discrepancy. 
However, these outliers correspond to isolated regions (the 
agreement between estimates is good in adjacent regions) 
and the overall x 2 = 44 for 36 d.o.f is satisfactory when 
these outliers are excluded. This suggests that the vi- 
gnetting correction is basically correct. We thus sum all 
spectra obtained for a given physical region and build up 
the temperature map presented in Fig. [|. 

There is no strong evidence of temperature variations, 
except for a cold area in the South-East (contiguous re- 
gions colder than average) and a hot area in the South- 
West. It is instructive to compare these temperature fea- 
tures with the X-ray image substructures (Fig. ||). The 
S/E cold region in the temperature map generally coin- 
cides with the filamentary substructure originating near 




Fig. 6. Excess emission over a j3 model overlayed on the tem- 
perature map. Isocontours are as in Fig. |l|. 
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Fig. 7. Temperature of each region of the temperature map 
with 90% confidence error bars. The temperature is plotted 
versus the distance to the brightest central galaxy NGC 4874 
and the data are splitted in three subsamples. Bottom panel: 
the S-E sector encompassing the filamentary structure towards 
NGC 4911 (open squares) and the S-W sector along the direc- 
tion towards NGC4839 group (filled squares). Top panel: the 
rest of the regions (filled circles). The angle (counterclockwise 
from west) defining the sectors are indicated in the figure. 

NGC 4911 and NGC 4921. It includ es the cold regions 
put into evidence by Donnelly et al. (1999) in that area 
(region 1 and part of region 20 in their Fig. 2), that we 
thus confirm. The hot regions in the S /W appear as a hot 
front perpendicular to the direction connecting the clus- 
ter center to the NGC4839 group, just ahead of the excess 
emission in that direction. This excess emission extends 
somewhat further to the North, where no specific temper- 
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Fig. 8. Radial temperature (top panel) and abundance (bot- 
tom panel) profiles. The rings are centered on NGC 4874. 

ature feature is apparent. However, the temperature map 
is specially noisy in that direction. The statistical signifi- 
cance of the temperature variations can be seen in Fig. [7] 
where we plotted the temperature of the various boxes ver- 
sus their distance to NGC 4874. We split the data in three 
subsamples: i) one S-E sector encompassing the filamen- 
tary structure towards NGC 4911 ii) one S-W sector along 
the direction towards NGC4839 group hi) the rest of the 
sample. The hot front (kT ~ 11 kcV) in the S-W located 
at about 15' from the center clearly stands out, as well 
as the colder region (kT ~ 6 keV) beyond 10' in the S-E 
sector. Otherwise the temperature does not deviate signif- 
icantly from the 8.25 ± 0.75 keV temperature range (less 
than ±9% variation). In particular, we see no evidence of 
the hot spot seen by Donnelly et al. (1999) 3' north of the 
NGC 4874 galaxy. We further extracted the spectrum (us- 
ing the central pointing) corresponding to thi s hot ASCA 
region: from Fig 1 of Donnelly et al. fll999|) we consid- 
ered a rectangular region of size 7.4' x 3.2' centered on 
a = 12 h 59 m 40 s , 5 = 28°00'30". We get kT = 8.4±0.4 keV 
consistent with the mean value and inconsistent with the 
ASCA value of Box 11 (kT = 12.7t 3 2 e keV). 

3.4. Temperature profile 

We further extracted spectra in concentric rings centered 
on NGC 4874 (all the data are summed). The region 
around NGC 4889 (40" in radius) and bright point sources 
were excluded. The temperature and abundance profiles 
are shown on Fig. ||. 

4. Discussion and Conclusion 

The dynamical state of the core of Coma has been much 
debated, in particular the nature of the merging unit(s) 
and their link with the dominant cluster galaxies. If there 
is a consensus that a merging group is associated with 
NGC4889, the situation of NGC4874 is less clear (see in 



particular Colless & Dunn 1996; Donnelly et al. 1999). 
As already proposed by Colless & Dunn (199E), our data 
suggest that NGC4874 is simply the central galaxy of the 
main Coma cluster, rather than being associated with a 
second subgroup in an early merging stage with the main 
cluster (as proposed by Donnelly et al. 1999). Several facts 
support this picture. First the remarkably homogeneous 
temperature distribution within the central 9 < 10' re- 
gion suggests the gas in that region is basically in a relaxed 
state. Second, there is no obvious evidence of a third peak 
in between NGC4889 and NGC4874, which would be asso- 
ciated with the 'true' cluster center. Actually, apart from 
the excess around NGC4889, the X-ray morphology can 
be classified as an offset-center cluster morphology (varia- 
tion of isophote centroid with scale) . The excess (size ~ 3' 
in radius) around NGC4874, when subtracting a (3 model 
(representative of the large scale morphology) is a natural 
consequence for this type of morphology. Moreover, part of 
the excess is certainly due to the contribution of the halo 
of the galaxy itself. The very significant drop in tempera- 
ture within V = 40kpc of NGC4874 (Fig. |) is natural in 
that context, as well as the increase in abundance, which 
could be due to an enriched ISM. 

It might be surprising that substructures survive in 
the gas density distribution (excess around NGC4889 and 
centroid shift for the main cluster) while the temperature 
distribution appears homogeneous. We must first empha- 
size that the spatial resolution and accuracy of the tem- 
perature and imaging data is not comparable. However, 
our results may also indicate that the gas simply follows 
the dark matter distribution. Dark matter substructures 
can survive for a long time after mergers , as in dicated by 
high resolution simulations (Moore et al. 1998 ). High res- 
olution hydrodynamic simulation, and more sophisticated 
morphology analysis, are essential to better understand 
this issue. 

At larger scale we do see evidence of recent merger 
activity. The cold filamentary structure in the South- 
East can be naturally explained by a merging group (see 
Vikhlinin et al. |1997| and Donnelly et al. |1999[). The posi- 



tion and extent of the cold substructure and the core prop- 
erties outlined above suggest that the merging group is as- 
sociated with NGC4911 and NGC4921, rather than being 
due to gas stripped from a group centered on NGC4874 
as proposed by Donnelly et al. ( 1999 ). Note that an ex- 
cess in the galaxy distribution is also observed around 
NGC4921/NGC4911 (Mellier et al. |1988[ ). 

Our analysis revealed for the first time a hot front in 
the South- West, just ahead of the excess emission that 
we see at the edge of the MOS mosaic and which extends 
further away towards NGC4839 (see Briel et al. |200lj ). It 
is situated roughly at the boundary of the group associ- 
ated wit h this galaxy, as defined from the optical (Colless 
& Dunn 1996j Fig. 9) and is perpendicular to the direc- 
tion connecting the center of Coma and NGC4839. This 
temperature structure is likely to be due to adiabatic com- 
pression, caused by the infall of matter associated with the 
NGC4839 group. The feature we find is indeed very sim- 
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ilarto the one displayed in Fig.6c of Schindler & Miiller White, S.D.M., Briel, U.G., Henry, J.P. 1993, MNRAS, 261, 



(1993), where at this time of the merger event no accretion L8 
shock has yet formed. To definitively characterise the fea- 
ture we need to know the temperature structure at larger 
radii, to quantify the transition conditions. It has already 
been noted that the NGC4839 group is located along the 
large scale filament connecting Coma and A1367 (e.g West 



et al. 1995). It is commonly thought that clusters form 
preferentially through anisotropic accretion of sub-clusters 
along large scale filaments. Our finding supports this sce- 
nario. The merger activity in that direction, particularly 
interesting for our understanding of cl uster formation, is 
further discussed in Briel et al. and Neumann et 

al. ( |200lD . 



Except for the very center as discussed above, the 
abundance is constant and the temperature radial pro- 
file is very weakly decreasing with radius. The slight drop 
beyond 10' is likely to be due to the cold S/W structure. 
The implications of this profile for the distribution of dark 
matter in the core will be studied in a forthcoming paper. 
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